INTERTWINING CERTAIN FRACTIONAL DERIVATIVES 



PIERRE PATIE AND THOMAS SIMON 

Abstract. We obtain an intertwining relation between some Riemann-Liouville operators of 
order a G (1,2), connecting through a certain multiplicative identity in law the one-dimensional 
marginals of reflected completely asymmetric a— stable Levy processes. An alternative ap- 
proach based on recurrent extensions of positive self-similar Markov processes and exponential 
functionals of Levy processes is also discussed. 



1. Introduction and statement of the result 
Consider for every a G (1, 2) the following operators acting on functions from 1R + to 1R: 



dx J T(2-a) — ' w dx 2 J x T(2 - a) 

for every x > 0. The operator is known as the left-sided Riemann-Liouville (RL) derivative 

of index a — 1 and D°L as the right-sided RL derivative of index a. Recently, RL derivatives have 
appeared quite often in various domains of analysis and probability and we refer to Chapter 2 
in p3] for a detailed account on them, as well as on other fractional operators. From the 
analytical viewpoint, RL derivatives extend in a non-local fashion the derivatives of integer 
order - see (2.2.5) in |14j . To name but one classical example of the occurence of RL derivatives 
in probability recall that 

Dt = D 2 o U 2 ~ a 

where D is the usual derivative and U 2 ~ a the potential associated to the standard (2 — a)-stable 
subordinator - see e.g. Exercise (1.6) in [4]. Consider now the operator 

A" = D%~ 1 o D, 

noticing first that it differs from D o D"" 1 = - see (2.2.29) in [14] for this latter equality - 
where D" is the left-sided RL derivative of index a which is analogously defined by 

d 2 f x /(tXx-t) 1 -" 



D " 11X) ~ r(2-a) dL 

Indeed, an integration by parts shows that for every x > 

(i-i) Km = + f[r^) /(0) - 
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Another fractional operator related to A" is the so-called Caputo a-fractional derivative which 
is given by 

C D%f[x) = —I— f ' f"{u){x-u?- a du, 
r(2 - a) J 

see section 2.4 in |14j . A similar integration by parts shows namely that 
(1-2) A«/(x) = C D%f[x) + f ^I^/'(0). 

For any a £ (1, 2), let V a be the multiplicative kernel acting on functions from R + to R by 

V«/(x) = E [/(^ a )] , 

where V a is a positive random variable having the density 

(-sin7ra)r- 2 (l + t) 
Va ^ ' ~ ir(t 2a - 2t Q cosavr + 1) 

- it will be checked soon afterwards that v a is indeed a density function on R + . Setting C 2 for 
the set of twice continuously differentiable functions R + — > R such that /' and /" are bounded, 
consider finally the following domain 

V = {/ G C 2 such that /'(0) = and 3 7 > 2 - a / lim x 7 (|/(x)| + = 0}. 

x— >+oo 

Observe that if / € X>, then necessarily 3 7 > 2 — a / lim x ._> +00 x 7 |/'(x)| = as well, which can 
be checked from the decomposition 

f(x + l) = f(x) + f'(x) + r + (x + l-y)f"(y)dy. 

J X 

It is also easily seen that D°_ and V a are well-defined on D, and it will be proved in the next 
section that A° is well-defined on V a (T>). Our main result is the following intertwining relation 
between A? and Z?": 

Theorem. For any f £ D, one has 

AlV a f = V a D a _f. 

Intertwining relations between Markov processes have some history and we refer to [12] 
for a probabilistic account, as well as various examples and applications. See also [18J for a 
particular analytical study. Though expressed in analytical terms, our result falls within the 
Markovian framework. It was namely observed in pQ that A™ (resp. D°t) is the infinitesimal 
generator of the spectrally positive (resp. spectrally negative) a— stable Levy process reflected 
at its running supremum. It might be interesting to mention that although the underlying stable 
Levy processes are in classical duality, this is no more the case for the reflected processes, so 
that our result can be viewed as a kind of intertwisted duality relationship for the latter. 

The proof of the theorem hinges upon the well-known criterion given in the Proposition 3.2 
of [6] and an identity in law connecting the running suprema of completely asymmetric stable 
processes which was recently obtained in [28], involving the positive variable T a with density 

(-sin7ra)(l -K 1/q ) 
na(t 2 — 2t cos air + 1) 
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Notice in passing that our above function v a is the density of the variable T Q , hence it is a 
density function. Some particular attention is paid to the functional domain upon which the 
intertwining relation holds and the set T>, which we borrowed from pQ, appears to be a reasonable 
and not too small candidate. With the help of some Suprun-type formulas for the resolvent 
of spectrally one-sided Levy processes which had been derived in |23j, we can also identify the 
Fellerian domains of A" and D" . In theory, those Fellerian domains yield an optimal formulation 
for the theorem, although they do not seem very tractable. 

In Section 3 we discuss another approach which consists in interpreting the stable process 
reflected at its infimum as the unique self-similar recurrent extension leaving continuously of the 
stable process killed when it enters the negative half-line. This identification had been roughly 
explained in Example 3 of |25| and here we can also check analytically that these two Feller 
processes have the same infinitesimal generators. This gives another proof, looking somewhat 
more unified, of all the results contained in the Appendix in pQ. An independent proof of the 
identity in law between suprema of completely asymmetric stable processes, which is the key- 
argument for the theorem, is also proposed, involving some closed formulas for the exponential 
functionals of certain Levy processes which had been established in [2H 122] . Though overall a 
bit lenghtier, we believe that this second point of view provides some unity to our interweaving 
relationship, which appears to be coherent with several apparently disconnected identities. 

To conclude this introduction, we stress that the positive random variable Zp with density 

(— sinvra) shi7r/3 
vr(a - l)(t 2 - 2tcosavr + 1) ~ vr/3(t 2 + 2tcos/3vr + 1) 
where (3 = a—1 E (0, 1), which can be viewed as a cut off Cauchy variable, has already occured in 
several distinct areas of the literature, especially through its power transforms Yg = z\ . See for 



instance Formula XI.11.6 in [30J for connections with the /3-fractional power of linear operators, 
Theorem 1.1 and Theorem 1.2 in [15] for mixture representations of the Linnik and Mittag-Leffler 
distributions of index f3, and more general geometric stable distributions, or Exercise 4.21 (3) in 
[8] which shows that Yp has the same law as the independent quotient of two standard positive 
j3— stable laws. The above variables V a and T a are less classical than Zp but one may of course 
wonder if they are not particular instances of a family of positive variables connecting suprema 
of general stable processes in duality, or a broader class of fractional operators than the one we 
consider in the present article. We plan to tackle this question in some further research. 

2. Proof of the theorem 

Let (Z,¥) be a spectrally negative stable Levy process of index a £ (1,2) starting at 0, with 
Levy density 

|y|-(a+l) 



so that the Levy-Khintchine formula reads 
(2.1) E 



e \z t 



e tx ° 



for every t, A > 0. See e.g. Chapters VII & VIII in [2] for an account on completely asymmet- 
ric stable processes. Setting St = sup{Z s , s < t} for the associated running supremum and 
introducing 

T x = inf{t > 0, Z t = x} = mf{t >0,S t = x}, 
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recall that T\ = S 1 a is a standard positive (l/a)-stable law, viz. 



E 



e _ A i/- 



for every A > 0. Denoting by It = inf{Z s , s < t} the running infimum, consider the reflected 
processes 

X t = S t - Z t and X t = Z t - I t . 
Notice that if Z = —Z is the dual process and if Yt = St — Zt and Yt = Zt — It are the 
corresponding reflected processes, then Y = X and Y = X. It is a basic fact from fluctuation 
theory - see e.g. Proposition VI. 1 in [2] - that X and X are Feller processes and we will denote 
by F x resp. F x their laws starting from x > 0. The infinitesimal generators L and L of X and 
have recently been expressed in [1], in three different forms. Recalling that by definition 

(00 , Tf( , v WHM ff ,, r t x [f(x t )]-f(x) 

(2.2) Lj\x) = lim resp. Lj{x) = hm 

t^-o f t->o t 

for every continuous function / : M + — )• R such that the limit in the right-hand side of the first 
(resp. second) equality exists uniformly, let us denote by Dom L (resp. Dom L) the set of such 
functions. Choosing then Riemann Liouville's form in Proposition A.l of [lj, one has T> C Dom 
L and 

Lf(x) = Dtf(x), x > 

for every / € D. Besides, one has 

Lf(x) = A%f(x), x > 

for every / 6 C fc 2 such that /'(0) = 0. The next proposition shows that A" is well-defined on 
V a (2?) so that the statement of our theorem makes sense: 

Proposition 2.1. One has V a {T>) C Dom L. 

Proof. Since v a (t) is of order t~( 1+a \ one sees by dominated convergence that for any / € T>, 
the function V a f is continuously differentiable on M + with bounded derivative 

(V a f)'(x) = E[V a f'(xV a )} 

- whence in particular (V a /)'(0) = 0, and twice continuously differentiable on (0,+oo) with 
second derivative 

(V a f)"(x) = E[V*f"(xV a )], x > 0. 
On the other hand, the right-hand side in the above equality might not be bounded when x — > 
because E[V^] = +oo. An easy change of variable shows however that for every / G T>, the 
quantity 

poo 

(xTAl) / {t l ~ a A I) f"(xt)dt 



JO 

remains bounded on (0, +oo) for some 7 < 2— a, so that (x 7 Al)(V Q /)"(x) is bounded on (0, +00). 
Hence, we need to show that in Proposition 4 of [TJ, the global boundedness assumption on F" 
can be relaxed into (x 7 A l)F"(x) bounded on (0, +00) for some 7 G [0, 2 — a). A persual of the 
proof shows that it is indeed the case: the relaxed condition changes namely the right-hand side 
of the second inequality in (4.11) therein into 



CaEy 



[ (Y s 2 ~ a + Yt a - r )ds 
Jo 
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for some positive finite constant C a , where stands for the law of Y (with their notation for 
Y, which matches ours) starting at y > 0. Using self-similarity, the expectation is then bounded 
by 

OE,[Y^]+E,[Y^]) < C an E y [Yt a ] 
for some other positive finite constants C' a ,C a ~. Notice now by the Feller property and the 
identification Y = X that the law of Y under is that of the process {Zt + (y V (—It)), t > 0} 
under P. Since < 2 — a < 1 < a, one has e.g. from Proposition VIII. 4 in [2] and the comment 
thereafter 

E^ 2 -*] = E[St a ] + E[(yV(-/ 1 )) 2 -] < +oo, 

which was already used in pQ. All of this shows that (4.12) remains unchanged for F under the 
relaxed condition, which entails exactly as in PQ that the limit in the right-hand side of (12. 2\\ 
exists pointwise for any / G V a (T>). To finish the proof, notice first by the Feller property that 
the transition densities of X form a strongly continuous contraction semigroup on the Banach 
space of continuous functions M + — > M tending to at infinity, to which belongs V a (D). Hence, 
one can apply the whole semigroup theory recalled e.g. in Section 31 of [26J. In particular, a 
result of K. Ito - see Lemma 31.7 in [26J shows that the limit in the right-hand side of (|2.2p is 
actually uniform for any / E V a (£>), whence / G Dom L as desired. 

□ 

Let us now consider the Fellerian domains of L and L. Setting Co for the set of continuous 
functions R + —> R tending to at infinity, the Feller property states that both functions 

x i y E x [f(X t )} and x ^ E x [f(X t )} 

are in Cq whenever f E Co. The Fellerian domains of L and L, which we denote respectively 
by T>(L) and T>(L) are made out of those functions in Cq such that the limits in (|2.2p exist 
uniformly - see e.g. Definition VII (1.1) in [24]. It is well-known from semigroup theory - see 
e.g. Proposition VII (1.4) in [24j - that 

(2.3) V(L) = U q (C ) and V(L) = U q (C ) 

for every q > 0, where U q , U q are the resolvent operators defined by 

f'OO roc 

U q f(x) = / e- qt E x [f(X t )]dt and U q f(x) = / e~ qt E x [f (X t )}dt. 
Jo Jo 

It follows from the definitions that 

DomL n C = V(L) and DomL n C = T>(L). 

The next proposition gives a full description of T)(L) and T>(L) in terms of the functions F a (x) = 
E a (x a ) and its derivatives, where E a is the Mittag-Lefner function of index a which is defined 
by 



E a (x) = Y,f 
For every f E Cq, introduce the further notation 



- ,an + 1) 
n=o v 1 



DC 



A/ = / e-Vf(y)dy. 



o 
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Proposition 2.2. One has V{L) = {14 XfF a (x) — F' a * f (x), f G Co} and 

( rx roa \ 

V{L) = [x^e- X J o K(y)f(y)dy + J (e^F^y) -F' a {y-x))f{y)dy, f G C j . 

Proof. Let us start with T>(L) = i7 1 (Co)- The resolvent density of X = Y killed when entering 
the half-line (a, oo) has been computed for every a > in |23j, Theorem 1 (i). Notice first 
that the notations therein yield Z^\x) = F a (x) and W^\x) = F' a {x). Besides, it was proved 
in Theorem 1 of [28j that the function F a — F' a is completely monotone, so that in particular 
F a (x) — F' a (x) — > as x — > +oo which, together with Formula 18.1 (10) in [9j, entails 

W {1 \a-y) _ y 

777- — - — — > e y as a — > +00 

ZW(a) 

for every y > 0. By monotone convergence, letting a — > +00 shows that the density of U 1 is 

v}(x,y) = e~ y F a (x) - F' a (x - y)l {y < x} , 

whence 

/ ;>oo \ px 

U l f{x) = / e~yf{y)dy)F a {x) ~ / F' a {x - y)f(y)dy = X f F a (x) -F^*f(x) 



for every / G Co an d x >0, which is the required identification for T>(L). 

Before identifying T>(L) = t/ 1 (Co), let us check that the two integrals in the definition are 
well-defined for every / G Co an d x >0 : the convergence of the proper integral comes from the 
easily shown behaviour F'^{y) ~ (a — l)y a ~ 2 /T(2a) as y — > 0, and the existence of the improper 
one is proved in combining several times formulas (6) and (43) in |llj . which yield 

aF' a {z) = e z + 0(z-( 1+a >) and aF^(z) = e z + 0(z~ (2+a) ) 

as z — > +00. The latter asymptotics also entail, with the notations of [23], 

W^ia-x) F' a {a-x) 

rfy = " n \ -»• e x as a ^+00. 

WX ] (a) F a(a) 

Since moreover W^^O) = ^(0) = 0, one obtains from Theorem 1 (ii) in |23| and monotone 
convergence the following expression for the density of U% : 

u^x^) = e~ x F^y) - F' a {y - x)l {y > x} 
whence, as above, the desired expression for V{L). 

□ 

Remarks 2.3. (a) If g = U 1 / G T>(L\ then g(0) = Aj and g is continuously differentiable on 
M + with derivative g'(x) = XfF' a {x) — F" * f(x), so that in particular g'(0) = 0. On the other 
hand, choosing / non differentiable shows that T>(L) contains functions which are not C 2 , as 
of course might be expected from the expression of L. See Chapter 2 in [TJ] for more material 
concerning the domains of Riemann-Liouville derivatives. 

(b) If g = U l f G T>(L), then 

/•oo 

9(0) = / {F^y)-F' a {y))f(y)dy 
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and g is continuously differentiable on R + with derivative 



■ : F'^y))f{y)dy 



so that in particular g'(0) = 0. Again, T>(L) contains functions which are not C 2 . 

End of the proof. We will use the inclusion (2) (3) in the Proposition 3.2 of [6]. We already 
know that X and X are Feller processes and it follows easily from the (1/a)— self-similarity of 
Z that they are also (1/a)— self-similar Markov processes. More precisely, one has 



{X b 



(X,¥ bx ) and (X b ,) 



for every b > 0, with X\ = bX h -a t and an analogous notation for X b . Since a / 1, this shows 
that they are not semi-stable Markov processes viz. 1-self-similar Markov processes - see (l.b) 
in [6]. However, a perusal of the proof of Proposition 3.2 in [6] shows that its statement remains 
unchanged when considering (1/a)— self-similar Markov processes for every a > and not just 
a = 1. 

We next show that the distribution of X is determinate under Po, with the notation of [6]. 
From e.g. Proposition VI. 3 in [2], the law of Xt under Po is that of St under P. Hence, taking 
the Laplace transforms, we need to show that if /, g € Cq are such that 



(2.4) 



E 



e~ xt f(S t )dt 



= E 


11 







-Xt 



g(S t )dt 



for every A > 0, then / = g. The latter is a basic property of (1/a)— stable subordinators but we 
will give some details for the reader's comfort. Recalling the notation T x = inf{t > 0, Zt = x} 
for every x > 0, we know that {T x , x > 0} is a standard (l/a)-stable subordinator since Z has 
no positive jumps. Besides, one has Leb (M + — U x >o(T x _, T x )) = a.s. - see e.g. the beginning 
of Section III. 5 in [2j, so that (f2T4l) entails 



E 



-XT X 



(1 



-\AT X \ 



.x>0 



E 



E 

,x>0 



g(x)e 



-xt x 



(1 



-XAT X 



The so-called Master's Formula - see e.g. Proposition XII (1.10) in [24] - yields then 



E ( 
' /o 

whence 



f{x)e 



-xt x 



-Xu 



,a+l 



j dx 


= E 











g(x)e 



-at, 



1 - e 



-Xu 



,a+l 



du dx 



f(x)e-^'~ x dx = I g(x)e- xl/ax dx 



—X x / a x 

^0 

and the required identification / = g by inversion of the Laplace transforms. 

Last, we see from Proposition VI. 3 in [2j and Formula (9) in [28] (beware the inverse notations) 
that 

X 1 = = T- 1 ^ x S 1 = V a x St = V a x Xt, 

where the identity T a = V a was mentioned in the introduction and follows from a change 
of variable. Putting everything together, Proposition 3.2 in [6j shows that A"V Q / = V a D"f 

for every / £ 'D(L) such that V a f £ "D(L). Supposing now that / £ V C DomI n Cq = 'D(L), 
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Proposition 12 , 1 1 shows that V a f £ Dom L and it follows immediately from dominated convergence 
that V a f £ Co- Hence, V a f £ DomL n Cq = T>(L) and we have shown 

&+V a f = V a D«f 

for every / £ V, as required. 

Remarks 2.4. (a) The interweaving relationship 

c D«V a f = V a D«f 

also holds over V because of the identity (|1.2p . Actually, since /'(0) = for every / £ 'D(L), 
the operator C D® coincides with the generator L of X. 

(b) As mentioned during the proof, Proposition 3.2 in [6] shows that 

A a + V a f = V a Dlf 

holds for every / £ 'D(L) such that V a f £ D(L). The domains T>(L) and T>{L) have been 
identified in Proposition 12.21 and allow rougher functions than T> but their formulations are 
unfortunately not very tractable, contrary to D. 

We close this section with an interesting complete monotonicity property for the Mittag- 
Leffler function E a . The latter is actually a direct consequence of Proposition 2 in [23J, but 
we present here a proof based on generators which is perhaps more transparent. Recall that a 
smooth function / :]0,+oo[— > R + is said to be completely monotone (CM) if 

d n f 

(-l) n ^L > 0, n> 1. 

By Bernstein's theorem, when f(x) — > 1 as x — > this is equivalent to the fact that / is the 
Laplace transform of a probability measure over R + . A classical result by Pollard states that for 
any a £ (0, 1] the function E a {—x) is CM, a property which does not hold anymore if a > 1 - see 
[27] for a proof of these latter facts and more on this topic. Like 1/x, the function x h-» E a (l/x) 
is CM by positivity of the coefficient in the series expansion - see Section 4 in |17j for further 
properties as well as general references on complete monotonicity. The following proposition 
shows a related property in the case a £ [1, 2]: 

Proposition 2.5. For any a £ [1,2] the function x t— >■ 1/E a (x) is CM. 

Proof. The case a = 1 is straightforward because E\{x) = e x and the case a = 2 follows from 
Levy's formula: one has 

= = E[e~ XT ] 

Ei (x) cosh yfx 

with r = inf{t > 0, \Bf\ > 2 -1 / 4 } and B a standard Brownian motion - see e.g. Exercise II 
(3.10) in [24J, so that \/E2{x) is CM by Bernstein's theorem. For the remaining case a £ (1,2) 
we will follow roughly the same arguments as Theorem 2.1 and 2.6 in |20j . Setting 

E%{x) = E a (qx a )f q {x) 

for every q > 0, where f q is some smooth function such that f q (x) = 1 for all x < 1 and f q (x) = 
for all x > 2, say, we see that E% £ 2?(A"). On the other hand, it follows from Lemma 2.23 in 
that 

A a + El(x) = qEUx) 
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for every x < 1. From Definition VII (1.8) and Exercise VII (1-24) in |24j - this shows that 
1 1 y e~ qT i At Ea(X T + At ) is a martingale under Pq, where 

T+ = inf{s > 0, X s > 1} = inf{s > 0, X s = 1} 

(recall that X has no positive jumps for the second equality). The optional sampling theorem 
entails 

1 = E [e-^ M El(X Ti+M )} -> E [e-^ + ^(X Ti+ )] = E a (q%[e~ qT ^} 

as t — > oo, by dominated convergence because X T + M € [0, 1] a.s. This completes the proof by 
Bernstein's theorem. 

□ 

Remark 2.6. It is plain from self-similarity, the Markov property and the absence of positive 
jumps that the variable is self-decomposable. By a well-known argument - see the proof of 
Proposition 2.6 in [20] for details - this shows that the mapping 

x i->- exp — [xE a {x)/E a {x)\ 

is also CM as the Laplace transform of a positive infinitely divisible random variable. 

3. Another approach with recurrent extensions 

The purpose of this section is two- fold. First, we will give another proof of the identifications 
L = D° and L = A?, viewing the reflected process X (resp. X) as a recurrent extension of the 
process Z (resp. Z) killed when entering the negative half-line, and using the classical expression 
of the infinitesimal generator of the unkilled stable Levy process. Actually the approach works 
for every strictly stable process Z such that \Z\ is not a subordinator, making it possible to 
retrieve the whole Appendix of pp. Second, we will derive a proof of the identity 

(3.1) X l = V a x X x , 

which is independent of |28j and relies upon closed expressions for the densities of exponential 
functional of certain spectrally negative Levy processes that had been carried out in |21|, [22] . 

3.1. Retrieving the generators of reflected stable processes. Let (Z,Q X ) be a strictly 
stable Levy process of index a S (0,2) such that \Z\ is not a subordinator, starting from 
We refer e.g. to Chapter VIII in [2J for details. The density of the Levy measure is 

v(y) = c + y _a-1 l {2/>0} + c_|y|" Q " 1 l {j/<0} , 

where c + ,c_ are nonnegative constants such that c + + c_ > 0. When a = 1 we suppose that 
Z is a symmetric Cauchy process viz. c+ = c_ = c > 0. Again, we will use the notations 
St = sup{Z s , s < t},It = ini{Z s , s < t},X t = St — Z t and X t = Z t — It- Setting also T = 
inf{i > 0, Zt < 0}, consider now the killed process 

R t = Z t l{T>t}- 

Again, since (T, Q x ) = (x a T, Qi) for every x > 0, one sees that (R,Q X ) is a positive (1/a)— self- 
similar Markov processes, viz. a Feller process taking values in R + (here, with as an absorbing 
state) and fulfilling the scaling property 

(R b ,Q x ) = (R,Q bx ) 
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for every b > 0, with the notation R\ = bRb-a t . As noticed in Example 3 of [25J, the reflected 
process (X,¥ x ) (with P x . defined analogously as in the preceding section) can be viewed as the 
unique self-similar recurrent extension of (J?, Q x ) leaving continuously. Roughly speaking, 
for every x > those two processes have the same law until the a.s. finite time T resp. 
T = inf{i > 0, Xt < 0} but is a regular boundary point for X, which is left instantaneously 
and continuously. The Feller process X has also infinite lifetime. See [10} 125] for precise accounts 
on recurrent extensions. 

We now identify the generator 1Z of the reflected process X viewed as a recurrent extension 
of R, retrieving in a unified manner all the results contained in the Appendix of pQ. Beware 
that since we consider the process reflected at its infimum, our notation is reverse to that of pQ. 
As in Section 2, we set (U q ) q >o for the resolvent of X. For every a £ (0, 1) and / G V, we use 
the same notations 

D«_f(x) = r fm ~ x)1 ~ a dt = -i r&t^n dt 

dx 2 J x T(2 — a) dx J x aF(— a) 

and 

V+nX) -dx 2 J r(2-a) dxj ar(-a) 

Proposition 3.1 (Bernyk-Dalang-Peskir) . For every f GT> and x > 0, one has 

C-/(0) 



i/a^l, and 



if a = 1. 



TZf(x) = r(-a)(c^D«f(x) + c + D«f(x)) + 



^/(•<-> '■( -^2 I f^ ] ^(r—^)'" + J{<)) 



Proof. Fix / G V and x > 0. For every M > 0, define /m over R in setting /m(x) = f{x) 
for every x > 0, Jm(x) = /(0) for every — M < x < 0, and letting — > smoothly as 
s — )• — oo. Then f^ £ except possibly at zero where its left and right second derivatives 

are bounded, and /m( x ) — > as |x| — )• +oo. Besides, with an abuse of notation, one can write 
TZf(x) = 7Zf'M(x) for every M > 0. Introducing the resolvent of R 

/•oo 

U q g(x) = / e-^ E x [g(Z t )l {T>t} ] dt 
Jo 

for every q, x > and g : M — >■ R + measurable, Theorem 2 (i) in [25] yields 

(3.2) ^/m(x) = \J q f M (x) + E x [e- qT ]U q f M ^). 

Recall from semigroup theory - see e.g. Exercise VII (1.15) in [21] - that 

lim q 2 U q f M {x) - qf M {x) = Hf M {x) = Kf(x), 

q — >oo 

which altogether with the notation T± = inf{t > 0, Zt > 1} entails 

Kf(x) = lim(q 2 U q f M (x) - qf M (x) + q 2 U q f M (0)^ x [e- qT }) 

= lim (q 2 U q f M (x) - qf M (x) + x~ a f(0)qE [ e -^}), 

q— »oo 
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the second equality being a consequence of self-similarity and the easy fact lim^oo qU q fM(0) = 
/m(0) = /(0). Proposition VIII. 4 in [2], a standard self-similarity argument and the Tauberian 
theorem quoted in [2] p. 10 give 

lim qE [e- qTl ] = re 

q— S-oo 

for some possibly vanishing, explicit constant re. This yields 

Kf(x) = lim (q 2 U q f(x) - qf(x)) + ^ 

and it remains to identify the limit on the right-hand side. Decomposing and changing the 
variable, one obtains 



/"OO f'OG 

q 2 V q f M (x) - qf M (x) = q 2 e~^ x [f M (Z t )}dt - qf(x) - / ie"* 

Jo Jo V */? 

By the same discussion as above, the Markov property at time T and the a.s. right-continuity 
of 1 1 — y Z^ at zero one has, recalling /m{ x ) = /(0) for every —M < x < 0, 

5«9> - £ (M) < 

E a [/ M (Z f/g )l {T < fM ] re/(0) , 

< hm sup ^ — L - L - < + e(M) 

g^+oo t/q x a 

for every t > and the same constant re as above, with e(M) — >• as M — > oo. By Fatou's 
theorem, this entails 



^x{fM{Z t /q)^{T<t/q}]\ (// 



«/(0) 



.X 1 



< limsup r te - t (^iM{z, lq )H? sl , } ]\ it s sfga + 



One the other hand, again from the resolvent equation, 

POO f'OO 

q 2 / e- qt K x [f M (Z t )]dt - qf(x) = q 2 / e-^ x \f M {Z t )\dt - qf M (x) -> £/a/(s) 

as g — )• +oo, where £ is the infinitesimal generator of Z. Indeed, one has /m G C|(R) except 
possibly at zero where its left and right second derivatives are bounded and /m(x) — > as 
|x| — > +oo, so that /m £ T>(C), as can be ssen readily from the proof of Theorem 31.5 in |26j . 
Supposing first 1 < a < 2, one has from the Levy-Khintchine formula 

£fAi(x) = / Um{x + y)~ fAi(x) - yf' M (x))u(y)dy 
Jr 

\f(x-y)-f{x) + yf(x))^ ri dy + ^-(/(0) - /(x)) + 
+ f°°(f(x + y)- f(x) - yf'(x))^-dy + e(M). 

Jo y + 
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Letting M — > +00, putting everything together and using the change of variable mentioned in 
[1 1 involving the assumption / G T>, one obtains 



TZf(x) = Cf(x) = T{-a){c-D a + f{x) + c + D«f(x)) + 



c-/(0) 



are- 



as desired. The cases a = 1 and < a < 1 are analogous and left to the reader. 



□ 



Remarks 3.2. (a) The above constant k can be identified as c+/a, see Lemma 3.1 in [5] and the 
references therein. The value of this constant does not play any role here, but it is interesting 
to note that it is exactly the same as the one extracted from the Levy-Khintchine formula in 
the above proof. 

(b) As in the proof of Proposition 12. 1\ it is possible to relax the condition / G T>. For example 
when a G (1,2) the global boundedness condition on /" can be changed into (x 7 A l)f"(x) 
bounded on (0, +00) for some 7 < 2 — a, and when a G (0, 1) the global boundedness condition 
on /' can be changed into (x 7 A l)f'(x) bounded on (0, +00) for some 7 < 1 — a. This is readily 
seen from the Levy-Khintchine formula and the proof of Theorem 31.5 in |26j . 

(c) With recurrent extensions, it is also possible to give an alternative proof to Proposition 
12.21 Suppose as in Section 2 that a G (1,2) and that Z has no positive jumps. By (|3.2p . a 
function is in T>(L) iff it can be written 

UJix) = Ui/(a;) + Exle-TpJiO) 

for some / G Cq. By self-similarity and Formula (7) in [28] - see also the references therein, we 
find first 

E x [e~ T } = Eofe-^ 1 ] = F a (x) - F' a {x). 



The term Ui/(x) can be handled with Suprun's formula. Specifically, letting a 
Theorem 1 of [3] and using the discussion made after Theorem 2 therein, we obtain 



K(y)f(x-y)dy 



+00 m 



Ui/(x) = I J e~yf(y)dy) F> a (x) - 
= X f F' a {x) - F^f(x) 
with the notation of Proposition 12.21 Last, we compute 



Uif(0) 



E 



e-tfiZt - I t )dt 



= E 
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-t 



f(S t )dt 



A 



where the last equality follows from the discussion after (|2.4p . paying here attention to the 
normalizing constants. Putting everything together yields the expression for T>{L) given in 
Proposition 12.21 The formula for T>(L) follows the same way, letting a,x,y —> +00 with a — x 
and a — y constant in Theorem 1 of [3] and identifying 



E 



(K(y)-K(y))Hv)dy. 



We omit the details. 
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3.2. Second proof of the theorem. In this paragraph we obtain a new proof of the identity 
(|3.ip which does not depend on the results of [28] but on Mellin inversion. More precisely, we 
will show that 

(3.3) Eo[X(] = E[V°] x E [Xf] 

for every s £ (1 — a, a), which is plainly enough to get (|3,ip . We start with the fractional 
moments of the random variable V a , a computation that could have been made directly by the 
residue theorem but since most of the argument was already carried out in [29] for some other 
purposes, we take the opportunity to shorten the proof. 

Lemma 3.3. For every s E (1 — a, a), one has 

sin(7r/a) sin(7rs) 



E[V* 



asin(7rs/a) sin(7r(l — s)/a) 

Proof. By equation (3) in [28] we know that the function 

= (-sinTra^q + t) 
J«W vr(t 2 « -2t a cosvra + l) 

is a probability density over IR + . The fractional moments of the corresponding random variable 
Y a can be computed with the help of the beginning of the proof of Proposition 4 in [29] and a 
change of variable: one finds 

_ sin(7r/q) sin(Trs) 

a a sm(irs/a) sin(7r(l + s)/a) 

for every s S (—a, a — 1). Notice that making s = — 1 entails 



/■OO 

/ v a (t)dt = Efy- 1 ] = i, 

Jo 



which shows that v a is a probability density with an argument slightly different from the intro- 
duction. Finally, the fractional moments of V a are given by 



EK1 = E[Y C 



s _, sin(-7r/a) sin(7rs) 



a sin(-7rs /a) sin(-7r(l — s) /a) 
for every s G (1 — a, a). 



□ 



To compute the fractional moments of X\ under Pq, we will need more material on recurrent 
extensions and exponential functional of Levy processes. With the notations of Section 2, let 
Z = —Z be the dual process and F x its law starting at x > 0. Introducing the stopping time 
T = inf{s > 0, Z s < 0}, consider the positive (1/a)— self-similar Feller process 

R t = Z t l {f>t} . 

The well-known Lamperti transformation [16] shows that the process defined ft = log R Tt for 
every t < T, with the notation 

T t = mils > 0, / R~ a dr > t 



o 
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is a Levy process starting at logx. Its Levy-Khintchine exponent tp which is defined by 

E[e~ x ^] = e^ (A) 

for every A > - recall that £ has no negative jumps so that the above expectation is finite, 
has been computed in [5] in terms of a certain improper integral. The next lemma gives a more 
tractable formulation in terms of Gamma functions. 

Lemma 3.4. With the normalization of Section 2, one has ip{X) = T{\ + a)/T(X). 
Proof. By Corollary 1 in [5] and Theorem 2.4 of [19], one has 



V>(A) = *(iA) = ? 7^^+ / (e- A ^-l + A(e^-l)l {e ,< 2} ) 

(a - l)r(-a) 7 



r(-a)(e» - 1) 1+Q 

+oo / y -A_ 1 + A(y _ 1} 



r(_ a )(i_j,)i+a y " y ^r(-a)(i- y )<* 



T(A) (a-l)r(-a) i„ Vr(-a)(l-y) 

where we have used several changes of variable. The last integral can be computed with the 
help of Formula (2.3) in |19| : one gets 

j Q { r[Z^ ) d y = -r(l-a)(2-a) a = l/(a-l), 
yielding the desired formula for ip. □ 

Remark 3.5. Supposing now that Z has positive jumps with the notations of Paragraph 3.1, 
an analogous simplification of Corollary 1 in [5] with Theorem 2.4 of [19] shows the following 
general formula for the Levy-Khintchine exponent of the Levy process associated to R by the 
Lamperti transformation: 

V>(A) = r(-a)(c_r(A + a)/r(A) + c+r(l-A)/r(l-(a + A))) 

for every A 6 (—a, 1). Setting 6 = inf{A > : = 0}, a simple analysis shows then that 

9 = ap = q(1 — p) where p is the asymmetry coefficient of Z, which can also be checked in 
considering the invariant function of R - see again Example 3 in |25| . 

Setting now 

6- = inf{A > : ip(—X) = 0} = 1 < a, 

Theorem 2 in [25] shows the existence of a unique recurrent extension for R leaving continu- 
ously, whose resolvent {U~) q >o is characterized by the formula 

-l/a poo 

(3.4) [7-/(0) = g / /(y)^- 2 !^^-]^ 

r(i - i/a)E[r /a x ] Jo 

for any positive measurable function /, where 



(3.5) L_ = / e*» 



oc 
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and £ is a spectrally negative Levy process with Laplace exponent 

r(a(A + l) - 1) 



^_(A) = V(«A-#_) = V(aA-l) 



T(aA - 1) 

As mentioned in Example 3 of [25], this recurrent extension is the process X, and with the help 
of this identification we can now compute the fractional moments of X\ : 

Proposition 3.6. For any s G (1 — a, a), one has 

sin(7r/a) sin(7rs)r(s + 1) 



E [Xf] 



aT{s/a + 1) sin(7rs/a) sin(7r(l — s)/a) 



Proof. Let us first connect the moments of X\ under Po to those of I_. Introducing the positive 
measurable function p s (t) = t s over R + , for every s G (1 — a, a) one has 

s/a+l 



s/a+l poo 

E °[ X il = rd/ n + U / e- qttS/a MXf]dt 
T(s/a + 1) J 



s/a+l roo q s / a+1 U- Ps (0) 
1 e- qt E [X t s ]dt ' 



T{s/a + l)J UL tJ T{s/a + l) 

where the second equality comes from self-similarity and the third from the fact that the resolvent 
of X is (U~) q >Q. From (|3,4p and after some simplifications, this entails 

F(l + ( s - l)/a)E[/Z (1+{s - 1)/a) ] 



E [Xi 



r(l - l/a)T(s/a + l)E[I. 



l/a-T 



for any s £ (1 — a, a). On the other hand, from Theorem 2.1 and Formula (2.1) in [21] - with our 
notation which entails 7 = l/a therein, see [22] for details - we know that the density function 
/_ of i_ has the alternate series representation 



/-W " a 2-A 1) ^ r(a(n + l)) V ' t>0 ' 



where C a is a positive constant to be determined below. This representation of the density 
prevents from computing the fractional moments of /_ by direct integration. Instead, one can 
use a so-called Mellin-Barnes integral representation of /_, which is obtained simply after a 
contour integration along a big half-circle in the half-plane x > —1 - see e.g. Section 3.4 in [13] 
for details. For any c G (—1,0) one has 

/-(*) : 





/•c+ioo 


2vri , 


/ c— ioo 


C a 


/•c'+ioo 


2vri , 


1 c'— ioo 



r(s)r(s - l/a)r(-s + 1 + l/a) 



T(as - 1) 



t~ s ds 



after a change of variable and taking any c' G (l/a, 1 + l/a). The inversion formula for the 
Mellin transform - see e.g. formula (3.1.5) in [13] - entails then 

E[I'_] = Cg T{S + l)r(S + 1 " l / a ^~ S + 



r(a(s + l)-l) 
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for every s £ (—1+1/ a, 1/a). Notice in passing, though we shall not need this, that making s = 
allows also to compute C a = T(a — l)/(r(l — l/a)r(l/a)). After some simple transformations, 
we finally deduce that for any s E (1 — a, a), 

T(l + (s - l)/a)r((l - s)/a)T(-s/a) sin(7r/a) sin(vrs)r(s + 1) 



aT{l — l/a)r(l/a)r(— s) aT(s/a + 1) sin(7rs/a) sin(7r(l — s)/a) 

□ 



End of the proof. The property that the law of X\ under Po is that of 5i = T l 1 ' a under P 
and a well-known, aforementioned moment formula for T\ entails 

f rx-i - r(s + 1) 
Eo[Xl] ~ r (s /a + i) 

for every s > —1, so that (j3.3|) simply follows from Lemma 13.31 and Proposition 13.61 
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